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Abstract 

O 

We consider the quantum mechanics of a system consisting of two identical, 
£SJ ■ Planck-size Schwarzschild black holes revolving around their common center 

of mass. We find that even in a very highly-excited state such a system has 
very sharp, discrete energy eigenstates, and the system performs very rapid 
transitions from a one stationary state to another. For instance, when the 
£SJ ' system is in the 100th excited state, the life times of the energy eigenstates 

are of the order of 10~ 30 s, and the energies of gravitons released in transitions 
. between nearby states are of the order of 10 22 eV. 
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Some years ago Hawking and others introduced an interesting idea about a possible 
spontaneous creation of virtual black hole pairs. jl|-|5| This process would be analogous 
to the spontaneous creation of electron-positron pairs in quantum electrodynamics, and 
the members of the black hole pairs would presumably be Planck-size objects. As it is well 
known, electrons and positrons may sometimes, in very favourable conditions, form a system 
called positronium, in which an electron and a positron revolve around their common center 
of mass. The possibility of a spontaneous creation of black hole pairs, then, gives rise to some 
very interesting questions: Could the black hole pairs sometimes form systems, analogous 
to positronium, where two microscopic black holes revolve around each other? What are 
the possible quantum states of such systems? What happens when the system performs a 
transition from one quantum state to another? Could one observe these transitions? 

These are the questions which will be addressed in this paper. We shall consider the 
simplest possible case where the members of the pair are Schwarzschild black holes, both 
having the same mass M. In general, the problem of quantization of a system containing a 
black hole pair presents immense difficulties: Indeed, even the classical solution to Einstein's 
equations describing a spacetime containing two black holes is unknown. However, if the two 
microscopic black holes are sufficiently far away from each other, they can be considered, 
to a very good approximation, as point-like particles moving with non-relativistic speeds. 
Moreover, the gravitational interaction between the holes can be described by the good old 
Newtonian theory of gravitation. Hence, the quantum mechanics of the black hole pair is 
described by the non-relativistic Schrodinger equation: 

2 fx r 

where /i := \M is the reduced mass of the system, and E is its total energy, r is the distance 
between the black holes. 

Our first task is to check whether the approximations we just made are valid. As it is 
well-known from elementary quantum mechanics, in stationary states the energy eigenvalues 
calculated from Eq. ([!]) are 

F - 1GM2 (?) 

where n — 1, 2, 3, ... . In these states the expectation values of r are 

2fr 2 or l r 1(1 + 1)^ 
~ GM$ + 2 ' (3) 

where I — 0, 1, 2, (n — 1) is the angular momentum quantum number of the system. Our 
approximations are valid if 

where 



2GM 

Rs := — (5) 
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is the Schwarzschild radius of a black hole with mass M. Comparing Eqs. (Q) and (|5]) we 
find that we must have 

—/=M < M Pl , (6) 
v n 

where 

[he 

M Pl :=J— ^2.2 x 10- 8 kg (7) 

is the Planck mass. Hence, either n must be very big or M must be much smaller than the 
Planck mass. The average velocity on the orbit is 

(v ave ) n = {—f-c. (8) 
Mpi n 

As one can see, the black holes move with non-relativistic speeds if Eq. (||) holds. In what 
follows, we shall always assume that M is smaller than, or equal to, the Planck mass, and n 
is very big. In other words, we shall consider microscopic black hole pairs in highly-excited 
states. 

Consider next the transitions between highly-excited stationary states. It follows from 
Eq. (Q) that the energy released when the system performs a transition where n is reduced 
by one is 

h n 6 n A Mpi 

For instance, if n is around ten and M is the Planck mass, the energy released is 10 -3 
times the Planck energy, or 10 6 J, which is about the same as the energy needed when an 
automobile is accelerated from rest to the velocity of 100 km/h. As one can see, enormous 
energies could, in principle, be stored in systems containing microscopic black hole pairs. 

Since the holes are assumed to be uncharged, one may expect that the main reason 
for transitions between stationary states is the quantum fluctuation of the gravitational 
field between the holes: the quantum fluctuation of the gravitational field perturbes the 
stationary states, a transition occurs, and a graviton is emitted or absorbed spontaneously. 
Since the black holes are assumed to be in a highly-excited state, and therefore relatively 
far away from each other, we may use the linear field approximation when investigating the 
perturbative effects caused by the quantum fluctuations of the gravitational field. 

The Lagrangian of a point particle moving in a weak gravitational field h^ u is 



L = -M^ + ^) — -. (10) 

Under the assumptions that |/t M „| <C 1, and the particle moves very slowly we may write 
Eq. (|10p, in SI units: 

L w -Mc 2 + -M(5 jk - h jk )x j x k - Mch 0j x j - -Mc 2 h 00 , (11) 
2 2 
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where j, k = 1, 2, 3 and x- 7 := dx^ jdt. Dropping the term — Mc 2 , which is a mere constant, 
we may infer that, in the center of mass coordinates, the Lagrangian of the black hole pair 
can be written, in effect, as 



1 



j x k 



L = ~n[5 jk - -h jk (-,t) - -h jk (--,t)]x 3 i 

-/ic 2 [M^) + M-^)], (12) 

where f is the vector joining the hole 1 to the hole 2, and x^s are defined such that r = x^ej, 
where e^'s are orthonormal basis vectors. 

Now, it is easy to see that the terms proportional to h 00 represent the gravitational 
potential energy of the system. These terms should give, when the system is in a highly- 
excited state, the Newtonian potential energy between the holes. The terms proportional 
to hoj and hj k , in turn, are related to the quantum fluctuations of the gravitational field. 
If we assume that transitions from a one stationary state to another are associated with 
spontaneous emissions or absorptions of gravitons, the only remaining terms, in addition 
to the Newtonian potential energy, are the terms proportional to hj k . That is because 
the "scalar" and the "vector" gravitons proportional to h 00 and h j, respectively, can be 
gauged away, and the physical gravitons correspond to the (jfc)-components of the field h^ u . 
Therefore, the Lagrangian of our system interacting with spontaneously emitted or absorbed 
gravitons is 

1 rr 1, J s 1, / r GM 2 , . 

Hence we find, under the assumption that hj k is small, that the Hamiltonian operator of the 
black hole pair interacting with gravitational radiation is 

v GM 2 
2/i r 

where 

Hint := ^\ r -,t) + y k (- r -,t)]p J p k , (15) 

and pj is the canonical momentum conjugate to x J . In these equations, objects equipped 
with hats are operators replacing the corresponding classical quantities. 

The transition rate Rj_>/ from the initial state \i) to the final state |/) can be evaluated 
by means of the Golden Rule 0. Evaluation of the transition rate involves quantization of 
the linearized gravitational field in the radiation gauge where h®^ = for every fi — 0,1, 2, 3. 
This has been performed in the Appendix 0, and we find that if exactly one graviton with 
angular frequency to is emitted in a transition, then 

R l ^ f = ^^\(f\cos(h-f)e ab p a p b \z)\ 2 dn, (16) 
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where k is the wave vector of the graviton, e ab its polarization tensor, and dQ is the solid angle 
into which the graviton emerges. The polarization tensors e^ ab and e^ ab corresponding to 
the two physical polarizations of the graviton have been chosen such that 

e (A')a6 € W = 26 XX (17) 
for every A', A = 1, 2. In the lowest order approximation we can write: 

R w = ^^K/k a W>l 2 ^- (is) 

Transitions corresponding to this expression may be viewed as gravitational analogues of 
the El transitions in atomic physics. We shall therefore call them as Gl transitions. 

The next task is to integrate Eq. (|T6| ) over all the possible directions into which the 
gravitons may emerge. This has been done in Appendix [B|. If the initial state \i) and 
the final state |/) are taken to be solutions to the Schrodinger equation ([[]) of the black 
hole pair, and one takes into account the fact that gravitons have exactly two independent 
polarizations, one finds that the integrated transition rate from the state |i) to the state |/) 
in Gl transitions is 

From this expression one can obtain the Gl selection rules for the black hole pair. In 
other words, one finds the transitions for which the integrated Gl transition rate ( |iTj| ) is non- 
zero. Since gravitons may be viewed as spin two particles, one might expect that for allowed 
Gl transitions the angular momentum quantum number I as well as the corresponding 
azimuthal quantum number mi could change only by zero, or by plus or minus two. A 
detailed investigation, which has been performed in Appendix 0, shows that this is indeed 
the case: the Gl selection rules are: 

Al = 0, ±2 , (20a) 
Am f = 0, ±2 . (20b) 

It is straightforward, although very laborious, to calculate the transition rates in allowed 
Gl transitions from Eq. (|T^|). This arduous task, with explicit expressions for transition 
rates in different Gl transitions, has been performed in Appendix D. For instance, if the 
principal quantum number n is of the order of 100, then it turns out that the transition 
rates in transitions between nearby states are of the order of 10 28 1/s, if Al = 0, —2, and of 
the order of 10 24 1/s, if Al = +2, provided that we assume that the mass M is equal to the 
Planck mass Mp\. Transitions between nearby states, however, are not necessarily the most 
favourable ones. In Figs. 1, 2, and 3 we have plotted the transition rates as functions of 
k, the difference between the initial and the final values of the principal quantum number 
n. One finds, for instance, that when Al = Ami — ~ 2, and Zj = m^ = 50, then the most 
favourable transitions are those where k ~ 20, and if, for the same values of U and m^, 
Al = and Ami — — 2, then for the most favourable transitions k ~ 10. If Al = +2 and 
Ami — 0, then k = 3 for the most favourable transitions. 
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As we have seen, Planck-size black hole pairs perform extremely rapid transitions from 
one stationary state to another. However, the problem is, whether the transitions are 
propably too rapid so that one cannot meaningfully talk about discrete energy eigenstates 
at all. To answer this question, one must calculate the life time of an initial state with big n, 
when all the allowed Gl transitions are taken into account. The life time of an initial state, 
in turn, is the inverse of the sum of the transition rates of all the allowed Gl transitions. 
Using the formulas of Appendix D, and assuming that transitions from higher to lower states 
are dominant, one finds that the life time of an initial state with n of the order of 100, is of 
the order of 

r 100 ~10- 30 s, (21) 

if M = Mpi. This result settles, within the approximations made in this paper, the question 
about the existence of sharp energy eigenstates. According to Heisenberg's uncertainty 
principle the natural line width of a state with life time r n is 

6E n ~ - , (22) 

and therefore 

5E 100 ~ 10" 4 J . (23) 
However, the energy difference between nearby states is, according to Eq. (||): 

E 100 - E 9g « 2 x 10 3 J , (24) 



and therefore: 

SE 



100 in-7 



-E-100 — Egg 



1Q- J . (25) 



In other words, the uncertainty of the energy eigenvalues of the system is very much smaller 
than the energy difference between nearby states. Hence, the energy spectrum appears to 
be discrete, at least as far as one can trust in the approximations made in this paper. The 
transitions between discrete energy eigenstates are extremely rapid, and the energies released 
are enormous. 

In this paper we have investigated microscopic Schwarzschild black hole pairs revolving 
around their common center of mass. Considering the holes as point-like objects interacting 
with the Newtonian gravitational force, we quantized the system and studied the stationary 
energy levels when the system is in a highly-excited state. We then calculated, by means 
of perturbative methods, the transition rates and life times in a certain class of transitions. 
These calculations were based on the quantization of the linearized gravitational field in the 
radiation gauge. We obtained, in the lowest order approximation, explicit expressions for the 
transition rates and evaluated the transition rates numerically when the pricipal quantum 
number n ~ 100, and the mass of a microscopic black hole is assumed to be one Planck mass. 
We found that the transition rates are of the order of 10 24_28 1/s for allowed transitions, 
and the life times of energy eigenstates are of the order of 10~ 30 s. Gravitons with energies 
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of the order of 10 22 eV, or even greater, are emitted in these very rapid transitions between 
discrete energy eigenstates. 

No doubt, one might be justified to have some feelings of suspicion towards the validity 
of the approximations on which this paper is based. Indeed, it may well be somewhat daring 
to apply Newton's ancient theory of gravitation to microscopic black holes! However, if one 
accepts the view that, in the classical level, Einstein's general theory of relativity is the 
correct theory of gravitation (so far we have no experimental evidence suggesting that the 
things could be otherwise), having Newton's theory as its non-relativistic limit, one is also 
forced to accept the "Newtonian approximation" made in this paper: Holes with M = Mpi 
move, if n ~ 100, with velocities which are of the order of 0.01c and the expectation value 
of their mutual distance is about 10 4 times their Schwarzschild radius. Of course, one 
could calculate general relativistic corrections to energy levels and transition rates but as 
far as one is interested in mere order-of-magnitude estimates, Newton's theory should be 
sufficiently accurate. Quite another matter is, whether the first order perturbative quantum 
theory of linearized general relativity is sufficiently accurate an approximation of quantum 
gravity for the evaluation of transition rates and life times. As far as the energies of the 
emitted or absorbed gravitons are well below the Planck energy, which is the case at least 
for transitions between nearby states when n ~ 100, however, one might be inclined to rely 
on the approximations made in this paper. 

One of the basic lessons one can learn from this paper is that enormous energies could 
be relased by means of the quantum effects of the gravitational field: The energies released 
in transitions between the energy eigenstates of the microscopic black hole pair are, even in 
the transitions between nearby states when n ~ 100, about fourteen orders of magnitude 
greater that the energies typically released in nuclear phenomena, and yet we are talking 
about a microscopic system. 

It appears to us that the only conceivable way how microscopic black hole pairs could 
be formed is a pairing of black hole remnants which are left as the end products of Hawking 
evaporation of primordial black holes. The existence of black hole remnants is a controversial 
question f7j, although it is possible to construct mathematically consistent quantum theories 
of black holes where the ground state energy of the hole is positive (see Refs. ||||). If such 
remnants exist, however, it is possible that some of them form pairs. Our calculations 
show that if the masses of the remnants are of the order of Planck mass, extremely energetic 
gravitons are emitted in very rapid transitions from one state to another. Some of these very 
energetic gravitons, in turn, may materialize into observable particles producing cosmic rays 
with very high energies. Measuring the amount of very energetic cosmic rays could therefore 
be used when trying to estimate the abundance of black hole remnants in the universe. If 
there are black hole pairs formed by black hole remnants with masses of the order of Planck 
mass, one expects to observe cosmic rays with energies of the order of 10 20 eV, or even 
greater. 

ACKNOWLEDGMENTS 

We are grateful to Matias Aunola, Markku Lehto and Jorma Louko for usefull discussions 
and constructive criticism during the preparation of this paper. 



7 



APPENDIX A: QUANTIZATION OF LINEARIZED GRAVITY IN THE 

RADIATION GAUGE 



The Lagrangian density of the linearized gravitational field h^ v can be written, when the 
Hilbert gauge condition 



^(/^-^r/i) = o (ad 



is used, as (c = 1) pi 



C = 7^(dxh^d x h» u - 2d fl h flu d x h l/X + 2& l h ilv d v h - d u hd u h) . (A2) 
In the radiation gauge we have, in addition 

h»° = V// = 0,1,2,3, (A3) 
and the Lagrangian density can be written in the form: 

C = -^(h mn h mn + dihnn&h™ - d n hd n h) , (A4) 

o47rG 

where h mn := dh mn /dt, and the Latin indecies take the values 1,2,3. The canonical momen- 
tum conjugate to h ab is therefore 

8C 1 

pfl6 _ OJ^_ = _J_ h ab (A5) 

dhab 32nG 

As the first step towards quantization we write the field h ab as a Fourier expansion: 

habit, r>EE <&W*> r-)al A) + ul(t, r)a^) . (A6) 

k A=l 

In this equation, the sum is taken over the wave vectors k and the polarizations A. e ab is 
the polarization tensor, gu and al are the Fourier coefficients, and the functions Uj: are 
orthonormal wave modes: 



u 



:(t,f) = Nte^'-^ , (A7) 



where Nt is a normalization constant, and u is the angular frequency of the graviton. 
Introducing periodic boundary conditions in a box with edge L: 

h ab (0,y,z,t) = h ab (L,y,z,t) , (A8a) 
h ab (x, 0,z,t) = h ab (x, L,z,t), (A8b) 
hab(x,y,0,t) = h ab (x,y,L,t) , (A8c) 

we find that the possible values of k x , k y and k z are 
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L ' 
2n 



k y — n y 



k 7 = n 



2vr 



(A9a) 
(A9b) 
(A9c) 



where n x , n y and n z are integers. Introducing, moreover, the inner product between wave 
modes, 



{ut,\uj:) := -i / dx dy dz {ut,u*r - UpU*r) 



"k lU, k "'k'"'k' ' 

together with the requirement that the wave modes are orthonormal: 

( u k'\ u k) = 

implies 

u % (t,r) = {2L a uj)- 



3, ,\-l/2 e i(k-r-ut) 



(A10) 



(All) 



(A12) 



Consider now the polarizations of gravitons. The Hilbert gauge condition (|A1|) , together 
with the radiation gauge condition (|A3|) implies 



f W k b 







(A13) 



for all a = 1,2,3 and A = 1,2. To satisfy this condition, we may choose the polarization 
tensors e^u such that 



e ab 6 



26 



XX 



(A14) 



for all A, A' = 1, 2. For instance, we may choose: 



ab 



(«S?) 




(Al5a) 



(A15b) 



In that case the gravitons propagate to z-direction. Using Eqs. ( | A 1 1|) and ( |A14j ) we can 



write the Fourier coefficients and a v r . Ajl as: 

k k 



,Wt 



,(A)t 



_ _ 1 (\)ab j dx j dy I dz{uj;h ab -u z h ( 

a w = * ( A)«6 r L dx r L dy f L dz ^*h ab -uZh ab ) 



(Al6a) 
(Al6b) 
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We now proceed to quantization. It follows from Eq.( |A5|) that the canonical equal time 
commutation relations between the field operators h ab and their conjugates p ab , 

[h ab (t,r'),p cd (t,f)} = th5 3 (?,r)5 c J d b , (A17a) 

[h ab (t^)Xd(t,r)} = [p ab (ty),p cd (t,f)} = 0, (A17b) 

can be written as: 

[h ab (t, 7), h (t, f)] = i32irGh6 3 {r', f)5 c a 5 d b , (A18a) 

xab ^cd 

[h ab {t,r'),h cd {t,r)] = [h (t,f>),h (t,f)]=0. (A18b) 

Using Eqs. ( |A14j ) and (|A16 ) we find that the commutation relations between the operators 
di and 5^ are: 



k k 



{&£>, 4 A)t ] = 16*GM y %- t , (A19a) 



a^,af] = [a<f*,aWt] =0 . (A19b) 

? creates and 

k k 

More precisely, they operate into the vacuum |0) such that 



Hence, a~ creates and annihilates a graviton with wave vector k and polarization A 



°k 



|0) = VlQnGh\l) , (A20a) 



4 A) |0> = 0, (A20b) 

where |1) is a one-graviton state. If we put everything here derived in together, we find that 
the operator H int of Eq. (|I5D takes the form: 

Hint = ^ E E(2L 3 c)" 1/2 e (A ^ cos(ifc • r){e^af + e^af^)p a p b . (A21) 

Z ^ l A=l 1 

Now, the Golden Rule || implies that if in a transition from the state |z) to the state |/) 
a graviton with energy flu, wave vector k and polarization A is emitted, the corresponding 
transition rate is 

R*-/ = ^M^L^- 11 ^ (A22) 

where 

p{Ef) '-=^ dn (A23) 
is the density of states close to the final state, and <El is the solid angle into which the 
graviton emerges. If the state \i) represents a zero-graviton state and the state |/) a one- 



graviton state, Eq.(A20) implies 



R, 



Guj 



or, in SI units: 



(f\ cos(lk-r)e^ ab p a p b \z)\ 2 cM (A24) 



27TC 5 ft/l 2 1 N 1 v 2 

which is the same as in Eq. (|T6|) 



^ ' (f\ cos(lk ■ r)e Wab p a p b \i}\ 2 dtl , (A25) 
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APPENDIX B: INTEGRATION OVER ALL DIRECTIONS OF THE 
GRAVITATIONAL RADIATION 



In this Appendix we shall calculate the integrated transition rate by integrating over all 
the possible diresctions into which the graviton may emerge. To evaluate the integral, we 
assume first that all the gravitons are the so called /i 12 -gravitons propagating to z-direction. 
Then the only non-zero components of the polarization tensor e ab are 

e 12 = e 21 = 1 . (Bl) 

Now, the idea is here that the components of the polarization tensor change when the 
direction of propagation of the emitting graviton changes. When the gravitons no more 
propagate to z-direction, we introduce a new Cartesian coordinate system K' which is ro- 
tated such that the gravitons always propagate, in frame K', along z-axis, and the original 
coordinate system K is kept fixed. The relationship between the orthonormal basis vectors 
e a and e' a of the frames K and K' is: 

e\ = sin e\ + cos e 2 , (B2a) 
e 2 = — cos 9 cos ei + cos 9 sin e 2 + sin 6 e 3 , (B2b) 
e' 3 = sin 9 cos t\ + sin 9 sin e 2 + cos 9 e 3 . (B2c) 

where 9 and are spherical angles. These relations imply that the polarization tensor 
corresponding to /i 12 -gravitons propagating to arbitrary direction is, in frame K , represented 
by a matrix: 

/ — cos 9 sin 20 cos #(2 sin 2 — 1) sin#sin0\ 
e— cos #(2 sin 2 — 1) cos 9 sin 20 sin 9 cos (B3) 
\ sin 9 sin sin 9 cos / 

To begin with we perform the summation in Eq. (|18|). Because the momentum operators 
commute with each other, we get: 



\(f\e ab p a p b \t)\ 2 dQ = J |e n (/|p 2 K) +e 22 W 2 |*> + 6 33 (/|p 2 K) + 2e 12 (f\p 1 p 2 \z) 

+ 2e 13 (f\pm\t) +2e 23 (f\p 2 p 3 \i)\ 2 dn (B4) 
After squaring and integrating the expression (|B4j) over all directions, we have: 



+ 8| W 2 p 3 |*>| 2 - (/&•>*</»> - OT 2 K>*</»> 

+ S(f\PiP3\i)*{f\P2P3\i) + 8(f\p 2 p 3 \ty(f\pm\t)) . (B5) 



In Appendix we shall calculate the transition amplitudes above in all details. If we now 
use the results (|D11|) , we get for the integrated transition rate a simple and nice expression, 
by taking into account that gravitons have two independent polarization states: 

^/=^i(/w>r, (B6) 

in the units where G = c = h = 1. 
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APPENDIX C: Gl SELECTION RULES 



In this Appendix we shall consider, in the lowest order approximation, the selection rules 
for the spontaneous emissions and absorptions of gravitons by the microscopic black hole 
pair. These selection rules we shall call Gl selection rules, in analogy to the El selection 
rules in atomic physics. The most convienient way of deriving these rules is to use the so 
called Wigner-Eckart theorem ||. That theorem concerns the irreducible tensor operators 
in the spherical basis. The Cartesian components pi t 2, 3 of the momentum operator can be 
written in terms of the standard components p-i, + i i of the irreducible momentum operator 
p of rank 1 as: 

Pi = —j= (p+i - P-i) , (Cla) 
P2 = (p+i + P-i) , (Clb) 

P3=P0- (Clc) 

To find the Gl selection rules, we must find all the allowed changes in the angular momentum 
of the microscopic black hole pair between the initial and final states. To know these changes 
one has to consider the transition amplitude given in Eq. (|B6|) . The transition amplitude 
can be written in terms of the standard components as: 

(f\fi\i) = U\\ (f+x -M-i -P-iP+i +P-i) \i) , (C2) 

where all the products between the standard components can be written, by using the 
definition of the tensor product, as: 

(f\fi\i) = l[(f\ £(1, !| J > M )W\Mi) - (/| EC 1 ' 1 ~ l \ J i M)[fp\jM\i) - 

1 J,M J,M 

(f\ 53(1,-1, 1, 1| J, M)\0\ JM \i) + (f\ 53(1,-1, 1,-1 1 J, M)\0\ JM \i)} , (C3) 

J,M J,M 

where [pp\jm denotes the JM-component of the irreducible tensor of rank 2, J corresponds 
to the eigenvalue of the total angular momentum operator J, and M is the corresponding 
eigenvalue of the projection operator J z . After taking the sum over J and M Eq. (|C3|) 
reduces to 

(M\i) = \{fMM - ^(f\moo\i) - ^(/|[^] 2 oN) + \{f\W\2 -2\i) , (C4) 



which, by Wigner-Eckart theorem, is equal to: 



2|„-\ _ I 1 \lf-rrif 



(MM = (-1) 



1 If 2 k 

2 I — rrif 2 mi 



(rif lf\\[pp\22\\ni k 



^(-X Om,)K '/HWoolk i,)- 
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Tff-X omjK U) + 



1 If 2 k 

2 I — Tflf —2 VTii 



n f h\\\Fi?\2 -2\\ni I, 



(C5) 



where and m/ are the eigenvalues of the z-component of the angular momentum operator 

at the final and at the initial states of the microscopic black hole pair, f ^ ft ^ | is a 

\mfMmiJ 

3j-symbol and (/||T^||i) is the reduced matrix element of the irreducible tensor operator of 
rank L. 

It is well known that the 3j-symbol vanishes unless the following two conditions hold for 
the eigenvalues of the angular momenta: 

\lf - L\ < k < l f +L (C6a) 
rrii - m f + M = . (C6b) 



These conditions, together with Eq. (|Cq), imply that the allowed transitions are the ones 
where Al = 0, ±1, ±2 and Ami = 0, ±2. However, there is a further constraint which 
comes from the parity conservation: Since the position operator x a is odd in reflections for 
every a = 1,2, 3, p a is also odd for all a. Therefore p a p b must behave as an even linear 
operator in reflections, and the states of the microscopic black hole pair must not change 
parity in Gl transitions. Since the parity of any state is given by (—1)', the value of / cannot 
be ±1. As the final Gl selection rules for the microscopic black hole pairs we find that the 
only allowed transitions in the lowest order approximation are the ones where 

A/ = 0, ±2, (C7a) 
Am; = 0, ±2. (C7b) 



APPENDIX D: TRANSITION RATES 

In this Appendix we shall derive explicit expressions for the transition rates in Gl transi- 
tions in the microscopic black hole pairs. In addition, we shall estimate the transition rates 
numerically in some well-chosen cases. Furthermore, we shall estimate the life time of the 
initial quantum state of the black hole pair. 

In a position representation the initial and the final quantum states, \i) and |/), re- 
spectively, are represented by the "hydrogenic wave functions", where the factor g 2 /47reo is 
replaced by GM 2 , where G is the gravitational constant and M is the mass of each hole | |11| : 

^i(^)WM), (Dl) 



GM 3 V (n-l-l)\ 
nh 2 ) 2n(n + l)\ 



exp 



-GM 
2nh 



-r) 

2 ' \nh 2 
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where 



L s k (x) 



E(-i: 



[s + k)\ 



iu M [s — m)\{k + m)\m\ 

is the associated Laguerre polynomial, and 



x 



(D2) 



y 

1 lm, 



-1 



(2Z + l)(Z-m z ) 



4tt(/ + m 2 ) 



; 1 [cos U ) exp{tmi( 



(D3) 



is the spherical harmonic function. P ; mi (x) is the associated Legendre polynomial. When 
calculating an expression for the transition rate, we will mostly be interested in the tran- 
sitions where the initial and the final states are in their maximal projection state where 
m\ i = U and m\ s = If. In that case we may use the result: 



p/(cos0) = (21 - l)!!sin'0. 



(D4) 



where the double factorial is defined as n\\ := n(n — 2)(n — 4) • • • • and 0!! = 1. 

In Appendix |C| we derived the Gl selection rules (|C7|). Because of these rules we must 
consider four different cases: We first calculate the three transition rates Ri^f where either 
Al = or Al = ±2 such that the initial and final states are maximal projection states, and 
then we consider the case where Al = and Ami — — 2. 



1. Al 



We first calculate the transition rate in a case where Al = — 2 = Ami ■ in other words we 
take the initial and final states to be \i) = \n, I, I) and |/) = \n — k, Z — 2, / — 2) (k = 1, 2, 3, . . .) 
such that I > 2 and n ^> 1. We can write these states in the position representation using 
Eq. (pi|) . To evaluate the integrals in Eq. ( P5|) , we must express the momentum operators 
p a in the spherical coordinates: 



Pi 

P2 
P3 



.*.{.„ ,9 cos 6 cos d 

-in smflcos</>— H — 

V or r oa 



,*_(.„. , d cos^sm^ d 

-in sm a sin <p— — I 

\ or r 



-ih cos 9 



d sin# d 



Or 




08 r sin 6 



r 89 



(D5a) 
(D5b) 
(D5c) 



After employing the expressions for the wave functions and for the momentum operators 
one gets a rather messy expression for the first transition amplitude in Eq. (|B5| ), which we 
still have to integrate (in the units where c = G = h = M = 1): 



(f\Pi\ i ) = (PifW) = {PifW) 

(2Z — 1)!!(2Z — 5)!! 
" SwnJ+^n - k) 1 



(21 + 1)(2Z - 3)(n - Z - l)!(n - k - Z + 1)! 
(n + Z)!(2Z)!(n-A; + Z-2)!(2Z-4)! 



1/2 
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nil n k poo ptt rlir 

y2 a m E b m i / drr 2 / d9sm9 / d(f) exp(z20) x 

m=0 m'=0 70 J ° J ° 

[sin 2 ' cos 2 (pR'^R'^r) + I sin 2 '" 2 cos 2 9 cos 2 <f>-Ri(r)R' f {r) - 
^sm 2l ~ 2 9sm2(i)R i {r)R'Ar) + sin 2 '" 4 cos 2 9 cos 2 <t>R'Ar)RAr) + 
/(j ~ 2) sin 2 '" 4 cos 4 9 cos 2 (j)RAr)R f (r) - ~ ^ sin 2 '" 2 sin 2(f>R'Ar)RAr) + 
sin 2 '" 4 sin 2 J R ? ,(r) J R / (r) 



2r 



where we have denoted: 
a m := (-l) r 

b m > := 
Ri(r) 



{n + l)\ 



(n-l-m - 1)!(2/ + m + l)!m!n r ' 
(n-Jfe + J-2)! 



(n - k - I - m' + l)\(2l + m> - 3)!m'!(n - k) r 
r 



2n 



R f (r) :=r /+m '- 2 exp 



2(n - k) 



(D6) 

(D7a) 
(D7b) 
(D7c) 
(D7d) 



and the prime ' denotes the partial derivative with respect to r. The 0-part of the integral 
is very easy to perform, and the radial part gives four different integrals which can be 
integrated separately by using the well-known result: 



POO 

/ dxx n exp(— ax) 
Jo 



ni 

,n+l 



(D8) 



Moreover, it is straightforward to show that the ^-integrals in Eq. (fD6|) can be reduced to 
the following integral: 



n n2l+ll]2 

dB sm 2l+1 9 ' ' 



o 



(2/ + 1)! 



(D9) 



When substituting all these integrals into Eq. ( P6|) , we get, for the transitions where Al = —2 
such that the system always remains in its maximal projection state: 



(M\i) 



1 



8n l+2 (n - k) 1 

n-l-l 

E 

m=0 
n-k-l+1 



21(21 -2){n + l)\{n -k + l- 2)!(n - I - l)!(n — k — I + 1)\ 



(21 + 1)(2/- 1) 



1/2 



X 



-1)' 



E 



(n - I - m - l)\(2l +m + l)\m\n r ' 

(-1)™' 



x 



7^ (n-k-l-rri + 1)1(2/ + m'- 3)!m'!(n - k) m ' 



x 



2n(n — fc) 
2n-k 



2i+m+m' 



(2Z +m + w! - 2)1(2/ + m + 2)m' 



2n-k 
2n{n — k) 
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,,/2/ + m'-4 2/ + m-l\ (21+m + m') 
(2l + m + m' - 1)! h 



2n 



2(n - fc) 



2(2n - fc) 



(D10) 



In a very similar manner one can show that the rest of the transition amplitudes in Eq. ( |B5| 
in this particular case, are given by: 



(f\pt\i) = -(f\m 

(Mi5 3 |i> = Vj = 1,2,3 
{f\Pm\i} = i(f\Pl\i) ■ 



(Dlla) 
(Dllb) 
(Dllc) 



The only difference between the rest of the amplitudes in Eq. ( |B5|) and the first ampli- 
tude ( |D10| ) comes from the 0-integration. 

The angular frequency u depends on the quantum number n and An = — fc as follows: 



2 [n 2 (n-k) 2 



(D12) 



As a final result, when the n-dependence of the angular frequency is taken into account, the 
transition rate for the microscopic black hole pairs becomes, in the case where Al = —2 = 
Ami: 



R 



1 



2/(2/ - 2)(n + Z)!(n -k + l- 2)!(n - / - l)!(n — k — l + l)\ 



24 \n 2 (n - fc) 2 / 

(n-l-l 



{21 + 1)(2Z - l)n 2l+ \n - k) 



21 



X 



E 



' (n-l-m - l)\(2l + m + l)\m\n r ' 



x 



n—k—l+l 

E 



7^ (n — k — l — m'+ 1)1(2/ +m' - 3)!m'!(n - fc) m ' 



2n(n — k) 
2n-k 



21+m+m 1 



(21 + m + m' - 2)!(2/ + m + 2)m' 



2n-k 
2n(n — fc) 



(2/ + m + m'- 1)! 



2/ + m! - 4 2/ + m - 1 \ (21 + m + m 1 



2n 



+ 



2(n - fc) 



+ 



2(2n - fc) 



(D13) 



If one prefers Si-units to the natural units, then the above transition rate should be multiplied 
by the factor Q grff L . 

As an example, let us now consider the binary black hole system consisting of two 
microscopic black holes with equal masses M = Mpi ~ 22/xg. In that case the transition 
rate between two neighbouring states from the state \i) = 1 100, 50, 50) to the state |/) = 
|99, 48, 48) is approximately: 



R 



fti 2.25 x 10 27 - 



and the corresponding life time of the initial state in this transition is 

1 



T 



R 



4.46 x 10 -28 



(D14) 



(D15) 



16 



When Am; 7^ —2, we may calculate transition rates from Eq. (|C5[) , since the reduced 
matrix elements do not depend on the value of the ^-component of the angular momentum 
operator. If Am; = +2 one can easily show that the transition rate is given by 

Rir+f=( ^ O - V f h l ^ ) 2 R^f> (D16) 

where -Rj-*/ is the transition rate related to the case Am; = —2 and the tilde ~ denotes the 
case Am; = +2. 

In the same way one can show that the case Am; = is related to the case Am; = — 2 
such that: 

Rir+f=( l ~ n - Vf I ^ ) '^/^ (D17) 

J \ —vdj \) rrii J V — W/ —2 m, / 

where -Rj^/ is the transition rate related to the case Am; = —2 and now the tilde "denotes 
the case Am; = 0. 



2. Al = 



Let us next consider the case where Al = such that the system always stays in its 
maximal projection state. In other words the initial state is \i) = \n, I, I) and the final state 
is \f) = \n — k, I, I). The position representation of these states can be read from Eq. (pip. 
In fact, the whole procedure to derive an expression for the transition rate in this case is 
analogous to the previous case where Al = —2. Therefore we just give the final result: 



R: 



1 

"3 \n 2 (n-k) 2 

(n-l-l 



/ + 1 \ (n + l)\(n - k + l)\ (n — I — l)\(n — k — I — 1)! 



E 



21 + 3 , 

_iyn 



21+4,, 



n — k) 



2Z+4 



X 



k 771 

n—k—l—1 



(n-l-m - 1)1(2/ + m + l)!m!n r ' 

(_lW 



x 



V 

^ (n-k-l-m' - 1)1(2/ + m' + l)!m'!(n - A;)" 1 ' 

2n(n-fe) 2Wm ' +21 
2n - k 



(2/ + m + m')! (/(m + m + I) + mm') 



6/ 3 



4/ 2 + 3/ + (m + m')(2/ 2 -4/)\ 2n - k 

2/ + 2 J 2n(n - k) + 



(2/ + m + m' + 1)! 



'/-m' + l / + m-l\ (2/ + m + m' + 2)!(n - fc) 



2n 



2(n - k) 



2n-k 



. (D18) 



For instance, if the initial and the final states are fixed such that the initial state is |i) = 
1 100, 50, 50) and the final state is |/) = 1 99, 50,50), then the transition rate between these 
states is 
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R 



« 1.44 x 10 28 - 



(D19) 



and the corresponding life time of the initial state in this transition is 

t « 6.94 x 1(T 29 s. (D20) 

When Ami 7^ 0, we cannot use Eq. (|C5|) to derive expressions for the transition rates; 
instead we have to perform another lengthy calculation. The difference between these 
cases comes from the definition of the associated Legendre polynomial, since we cannot use 
Eq. ( p4| ) , which holds only for the maximal projection states, in the cases where Ami ^ 0. 
Because our system is not necessarily in a maximal projection state we have to use the 
general definition for the associated Legendre polynomials: 



2H\ 



i=0 



iv 



3=0 



(D21) 



As this definition comes to effect in Eq. ( pi|) , the transition rate still is given by (|Bq). One 
can show, after performing some integrals, that when Al = and Ami — — 2, the transition 
rates of our system obey, for large n, the following expression: 



R 



1/1 1 \ (2/+ l) 2 2 2m 'mi! 2 (Z -mi)\(l -mi + 2)\{n + l)\ 

"24 \n? ~ (n - k) 2 ) n 2l+A (n - k) 2l+A (l + m{)\{l + mi - 2)! 

{n—l—l ^ — fc — Z — 1 I I 
E E EE a «An' x 
m=0 m'=0 s=0s'=0 

hK ss> T l ss , + I 3 (-K SS ,T 2 S , + K SS ,T*\ + I 2 (-K SS ,T* S , + K ss ,T^ s )j + 

y^ 2 

— „ , , _I_ If .T" 8 I 



X 



h I K SS ,T® S , - K SS ,T] S , + K SS ,T* S , 



where we have denoted 



b m ' 



h 
h 
h 
h 

Mr) 



(n - I - m - \)\(2l + m + l)\m\n m ' 

(-l) m/ 

(n-k-l-m 1 - l)\(2l + m' + l)!m'!(n - k) 7 

,2 D / 



drr 2 B^(r)R! f (r) . 
/ drrR'^Rfir) , 

« 

drrRi(r)R'f(r) , 
drRAr)Rf(r 



r l+m exp I — 



2n/ ' 



(D22) 

(D23a) 

(D23b) 

(D23c) 
(D23d) 
(D23e) 
(D23f) 
(D23g) 
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Rf (r). = r^e,^- w - ¥j j, (D23h) 



1 ss' 



/-i \2l— s— s' l+mi-1 l+mi-3 

K — ww^W^f. S l2s ' l{ fi |2s,_il ' (D23l) 

T 1 2'+«'-'( a + *'-/)! 

ss " (2s + 2s'-2Z + 3)(2s + 2s'-2Z + 1)! ' 1 JJ 

2 S+S "'(s + ^-O' ( 1 + 2 / _ _ A (m3k) 

(2s + 2s'-2Z + 1)! V2(2s + 2s'-2/ + 3) / 

ss " 2 s + s '-'+ 2m Ks + s'-/-l)!(2m i -l)m / ! 2 ' 1 ; 

7l, := ; Vt^ - + 2s' - mi - I + 2 , D23m 

2m t (2mi - 2)(2s + 2s' - 2l)\ \2mi - A / ' v ; 

_ (2-m f )(2m ; )!2^- 2m '( g + ^-Q! 

: " m^{2s + 2s'-2l + l)\ ' [UZ6n) 

t6 _ 2 a+a '- l (s + s'-l)\ f m-2 

2(2s + 2s' - 2l)\ W2s + 2s / -2Z + l)(2m i -2) 



2s' -I -mi + 2 (mi - 2)(2s - I - mi) 



(2s + 2s' -21 + 3) (2s + 2s' - 2/ + 1) mj(2m; - 2) 

(2s - I - mi)(2s' -l-m l + 2)\ 

m,(2s + 2s'-2/ + l) ) ' 
(2m / )!2 s+s '-^ 2mi (s + s' - l)\ f mi(2 - mi) 



+ 



(D23o) 



mi\ 2 \(2s + 2s' - 2l)\(2mi - 5)(2m, - 3)(2m z - 1) 

(2 - m{)(2s - I - mi) m t (2s' - I - m t + 2) ^ 



(s + s' - /+)!222s+2s'-2/ (2 m - 3)(2m, - l)(2s + 2s' - 2l)\ 

mi(2 - mi) ' 

(2mi - l)2 2s+2s '- 2l (s + s' - l)\ 2 (2s + 2s' -21 + 2) 



(D23p) 



mR (2-mi)2 s+s - l (s + s' , n s 

TL := ^ — ^ ^ — -L . (D23q) 

2(2s + 2s' -2Z + 1)! V HJ 

When Ami — +2, we may use Eq. (|U5|) to give an expression for the transition rate 



If 2 k V / If 2 U 



-2 



where -R^f corresponds to the case Al = and Am; = —2. 



3. Al = +2 

The only case we are left to consider is Al = +2. If, for now, we only allow the transitions 
where Ami = +2, then all the calculations, except for the angular part, proceed in a manner 
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very similar to the first case considered in this Appendix. However, the angular integrals 
are straightforward to evaluate. Therefore we just present the result of the calculation. The 
transition rate in the case An = —k, Al = +2 and Am; = +2 for large n is: 



(2l + 3)(2l + 5)n 2l + 4 (n 




(21 + 2)(2Z + 4)(n + Z)!(n - Jfe + I + 2)! 



-l) m 
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l)\(2l + m + l)\m\n m 

' 2 n( n -A;) 2/+m+m,+4 



3)!(2Z + m' + h)\m'\(n — k) 1 



In-k 



x 



(2Z + m + m' + 4)! (l + m' + 2 



2n 



An 
k 



k 



2n 



+ 



I + m 1 \ . , , 

+ - ] (21 + m + m +3)1 + 



(l + m! + 2)(l + m) 



2n-2k ' 2 / 
(2Z + m + m! + 2 



(D25) 



2n(n — k) 

As an example, let us now consider a transition between two neighbouring states from a 



state \i) 



|100,50,50) to the state |i) -- 

Ri^f w 8.30 x 10 23 
and the corresponding life time is 



99,52,52). We have: 
1 

s 



r w 1.20 x 10" 24 s. 

When Am; 7^ +2, we may calculate transition rates from Eq. 
can easily show that the transition rate is given by 

2 / , „ , \ -2 



(D26) 
(D27) 

If Am; = —2 one 



R; 



*f 



h 
-rh f 




If 2 li 
-rrif 2 rrii 



Ri 



and in the same way one finds that the case Am; 
Am; = —2 in the following manner: 

2 



is related to the case Al 



R; 



*f 




If 2 1^ 

—771/ 2 mj 



-2 



Ri- 



(D28) 
+2, and 

(D29) 



These relations conclude our considerations in these particular cases. 

To answer to the problem concerning the discreteness of the energy spectrum of the 
microscopic black hole pairs, we have calculated the life time r of the initial state of the 
system when all the possible transitions from a fixed initial state are taken into account. We 
have chosen the initial state to be |z) = 1 100, 50,50). This initial state yields a summation 
of about three hundred transition rates, when all the allowed Gl transitions are taken into 
account. A numerical estimate for each transition rate can be calculated by using the 
expressions for transition rates given in this Appendix. The life time of the initial state is 
obtained as an inverse of the sum of all the possible transition rates. We have found that 
when all the transitions are taken into account, the life time of the initial state is 



1.14 x 10 



-3CL 



(D30) 
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FIGURES 
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FIG. 1. The integrated transition rate Ri^f as a function of k, which is the difference between 
the initial value rii and the final value nj of the principal quantum number n. In this figure we 
have Al = Ami = —2. The initial values of n, I and m; are, respectively, rii = 100, k = 50 
and mi- = 50. In all of the figures Ri^f has been plotted in the units of x 10~ 15 , where 

f pl : = ~ 5.4 x 1Q- 44 s i s the Planck time. 



0.16 
0.14 
V 0.12 
5 0.1 
2 0.08 
| 0.06 ■ 
0.04 ■ 



5 10 15 20 25 30 35 40 45 50 



FIG. 2. R{->f as a function of fc, when = 100, ^ = = 50, Al = +2, and Am; = 0. 
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FIG. 3. Ri^f as a function of fe, when rii = 100, k = = 50, Al = 0, and Ami 
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